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Bacterial biofilms are usually assumed to originate from individual cells deposited on a surface. 
However, many biofilm-forming bacteria tend to aggregate in the planktonic phase so that it is 
possible that many natural and infectious biofilms originate wholly or partially from pre-formed cell 
aggregates. Here, we use agent-based computer simulations to investigate the role of pre-formed 
aggregates in biofilm development. Focusing on the initial shape the aggregate forms on the surface, 
we find that the degree of spreading of an aggregate on a surface can play an important role in 
determining its eventual fate during biofilm development. Specifically, initially spread aggregates 
perform better when competition with surrounding unaggregated bacterial cells is low, while initially 
rounded aggregates perform better when competition with surrounding unaggregated cells is high. 
These contrasting outcomes are governed by a trade-off between aggregate surface area and height. 
Our results provide new insight into biofilm formation and development, and reveal new factors that 
may be at play in the social evolution of biofilm communities. 


INTRODUCTION 

Surface-attached communities known as biofilms are 
believed to be the predominant mode of existence for bac¬ 
teria in many environmental settings [1]. Understanding 
how biofilms establish and grow is also clinically impor¬ 
tant given their ubiquity in medical implant infections [2], 
chronic wounds [3], and in the respiratory tracts of cystic 
fibrosis patients [4]. In the clinical context, biofilm com¬ 
munities often show enhanced virulence [5], resistance to 
antibiotics [6], and resistance to the host immune system 
[7]. These features may be associated with the spatial 
structure of the biofilm, which not only affects material 
transport, e.g., penetration of nutrients/antibiotics, but 
is also associated with differences in metabolism and gene 
expression among cells within the community [8, 9]. 

In the canonical picture of biofilm development, indi¬ 
vidual cells land on a surface, attach and proliferate to 
form first micro-colonies and later 3-dimensional struc¬ 
tures [10]. However, bacteria are also known to form 
dense aggregated clumps when they are grown in liquid 
(planktonic phase) [11-13]. Moreover, cells often disperse 
from existing biofilms as clumps of aggregated cells. Thus 
it is very likely that when a biofilm forms, some cells may 
arrive on the surface already in an aggregated state. In 
support of this view, evidence exists for the seeding of in¬ 
fections by pathogenic bacteria already in an aggregated 
state [14, 15], and bacterial aggregates are abundant in 
cystic fibrosis [4, 5] and tuberculosis [16] infections. 

Having arrived on the surface, e.g., a plant leaf [17], a 
surgical implant [2] or an industrial component [18], it is 


to be expected that cells within a bacterial aggregate will 
have to compete during biofilm development, both with 
other aggregates and with initially non-aggregated cells, 
to which they may or may not be genetically related. 

We take a first step towards understanding the role of 
pre-formed aggregates in biofilm development by investi¬ 
gating this competitive process, using agent-based simu¬ 
lations. Such simulations, in which the spatial structure 
of a biofilm emerges from local interactions between in¬ 
dividual cells, have become a staple tool for investigating 
biofilm structure and dynamics [19-21], as well as so¬ 
cial evolutionary aspects of biofilm development [22, 23]. 
Using this approach, we determine how a pre-existing 
aggregate of bacteria impacts the spatial structure of a 
biofilm, both in the presence and absence of competing 
unaggregated bacterial cells. 

Our main focus here is on the role of the initial shape 
of the aggregate. It is well known that bacterial inter¬ 
actions with a surface depend on features such as extra¬ 
cellular polymeric substances (EPS), presence of cell sur¬ 
face appendages (such as pili), and cell surface charge, 
which are species- and strain-dependent [24]. Moreover, 
soft-matter science has established that the nature of 
material-surface interactions can drastically affect the 
shape of fluid or semi-fluid droplets on surfaces [25]. It 
is therefore reasonable to suppose that in some circum¬ 
stances, bacterial aggregates will spread out in contact 
with a surface, while in other scenarios, aggregates will 
adopt a more compact configuration. Here we investi¬ 
gate the biological consequences of aggregate shape in 
the seeding of biofilm growth. 
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FIG. 1. Our simulation set-up. (a) Schematic representation 
of bacterial aggregates (green) which are initially spread on 
a surface to varying extents. The schematic also shows sur¬ 
rounding, competing, unaggregated cells (red). Top- Rounded 
aggregate, 0 = 180^; Middle- Semi-spread aggregate, 0 = 90°; 
Bottom- Spread aggregate with ^ ^ 0°. Note that the size 
of the aggregates (in terms of number of bacteria) is approxi¬ 
mately equal, (b) Growth of the spread {0 = 5°), semi-spread 
(0 — 90°), and rounded aggregate (0 — 180°) populations over 
the course of our simulations in the absence (p = 0 cell pm“^) 
and presence of competition (p = 0.5 cell pm”^). For clarity 
the error bars, representing the standard deviations, are only 
shown for the final data points. The standard deviations at 
these points are maximal. 


Simulating the development of biofilms initiated from 
initially spread or rounded aggregates, we find that the 
initial configuration of a bacterial aggregate on a sur¬ 
face is crucial in determining its eventual fate within 
the biofilm. In the absence of competitor cells on the 
surface, aggregates that maximise the extent to which 
they initially spread on the surface perform better than 
rounded ones because their initial access to nutrients (in 
the surrounding media) is greater. However when faced 
with strong competition from neighbouring unaggregated 
cells, initially rounded aggregates perform better at long 
times, despite the fact that the rounded aggregate shape 
has a smaller surface area and hence a reduced expo¬ 
sure to nutrients. Importantly, we show that in an ini¬ 
tially rounded aggregate, cells at the top of the aggregate 
proliferate at the expense of cells in the aggregate cen¬ 
tre. This has interesting possible consequences for social 
evolution given that cooperation within clumps of aggre¬ 


gated cells has been suggested to be a stepping stone in 
the evolution of multicellularity [26, 27] 

Our study highlights the effects of nutrient gradients 
and bacterial aggregate shape on long-term biofilm devel¬ 
opment. Our work reveals that these factors alone can 
produce a trade-off between nutrient access and compe¬ 
tition, with the balance between these factors depending 
sensitively on aggregate shape. While the link between 
biofilm spatial structure and nutrient access has been 
highlighted in many other studies [8, 23, 28-30], our work 
is the first to focus on the role of preformed aggregates 
in this context. Our study should help to decipher the 
role of pre-formed aggregates in biofilm infections. More 
generally, our findings emphasise the need to consider 
pre-formed aggregates in our current understanding of 
biofilm development. 


METHODS 

In this study, we use agent-based computer simula¬ 
tions to model the growth of a biofilm on a surface, start¬ 
ing from initial configurations of bacterial cells like those 
shown in Figure 1(a). In our simulations, an initial ag¬ 
gregate of cells (shown in green in Figure 1(a)), adopting 
a particular shape, seeds an inert surface, and may com¬ 
pete with surrounding unaggregated cells (red in Figure 
1(a)). Note that the red and green bacterial cells differ 
only in the manner in which they are initially arranged 
on the surface. At the start of our simulations, the “red” 
bacteria are distributed at random across those parts of 
the surface not occupied by the aggregate (see Sections 
SI to S5). To vary the extent of competition between the 
aggregated and unaggregated cells, we varied the initial 
cell density (number of cells per unit length of surface) of 
the unaggregated “red” cells (see Sections SI to S5). As 
a control, we also ran simulations in which the aggregate 
grew in the absence of the unaggregated cells. 

The focus of this work is on the shape of the initial 
cell aggregate. Figure 1(a) illustrates three different sce¬ 
narios, in which the cell aggregate adopts a compact 
rounded shape (top), spreads out on the surface (bot¬ 
tom), or adopts an intermediate shape (middle). In each 
scenario, the dimensions of the aggregate are adjusted so 
that the number of cells within the aggregate remains the 
same (see Sections SI and S2). 


Aggregate Shape Characteristation 

To characterise quantitatively the aggregate shapes in 
the different scenarios shown in Figure 1(a), we define 
the “aggregate-surface angle” 6>, which is the angle that 
the initial aggregate makes with the flat surface. A small 
value of ^ (^ ^ 0°) describes an initial aggregate config¬ 
uration that is spread on the surface, whereas a large 
value oi 0 {0 ^ 180°) describes an aggregate that is 
rounded. Given that the total number of cells in the 
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initial aggregate is fixed, 0 also encapsulates an interplay 
between the initial surface coverage of the aggregate and 
its initial height; with increasing 0^ the surface coverage 
decreases whereas height increases. In soft matter sci¬ 
ence, an analogous parameter is often used to describe 
wetting interactions between liquid droplets and surfaces 
[25]. A similar approach has recently been applied to the 
surface-spreading behaviour of eukaryotic cell aggregates 
[31]. From a phenotypic perspective, 0 is related to the 
nature of the interactions between cells in the aggregate 
and between cells and the surface, and thus could be tun¬ 
able by biological regulatory processes, or by evolution. 
We performed simulations for a range of 0 values between 
5° and 180° (Figure 1(a)). 

The aggregate configurations that we used to initiate 
our simulations were generated by “transplantation” of 
circular segments from simulation snapshots of pre-grown 
biofilms (see Sections SI and S2). This procedure proved 
preferable to other initialisation methods as it ensures no 
overlap between individual bacteria and enables the gen¬ 
eration of different aggregate shapes of the same number 
density (^100 cells per unit area). By varying the radius 
of the circular segment we are able to ensure that each 
aggregate contained ^100 cells that were initially spread 
on the surface to different extents. To ensure statistical 
accuracy, four different configurations were generated for 
each aggregate shape, defined by its value of 6>, and for 
each of these configurations, five simulations were per¬ 
formed using different random number seeds. Changing 
the random number seed affects the order in which in¬ 
dividual bacteria grow and divide, and also changes the 
locations of the unaggregated cells on the surface sur¬ 
rounding the aggregate. A total of twenty simulations 
were therefore performed for each value of enabling us 
to sample both variation in the configuration and the or¬ 
dering of cell updates (Section S5). Increasing the num¬ 
ber of repeated simulations did not affect our results. 

In common with many other biofilm simulation studies 
[19, 22, 23, 32], our simulations were performed in two 
dimensions for the purposes of computational efficiency. 
We have verified, however, that our key findings are re¬ 
produced when we use 3D simulations (see Sections S3 
and SIO). 


Simulation Implementation 

We use the agent-based microbial simulation package 
iDynoMiCs [33] to model biofilm growth, starting from 
configurations such as those shown in Figure 2. In these 
simulations, individual bacterial cells are represented as 
spherical agents, which grow and proliferate conditional 
on the local nutrient concentration, and “shove” each 
other apart to relieve local stresses within the biofilm. 
The order in which cells are selected to grow and divide 
is random during each global time-step of the simulation, 
as is the direction of cell division. In our simulations, the 
initial distribution of surrounding competitor cells on the 


surface is also random. The simulations use a spatial grid 
to track the local nutrient concentration field. Nutrient is 
assumed to diffuse towards the biofilm from above, with 
the concentration being fixed to a bulk value in a layer 
far from the biofilm. Within the biofilm itself, nutrient 
diffusion is hindered relative to the region outside the 
biofilm. Nutrient consumption by the bacterial cells leads 
to local gradients, which can have a strong impact on the 
structural features of the growing biofilm [8, 23, 28-30]. 
Periodic boundary conditions are imposed on both the 
nutrient concentration field and the particle coordinates 
in the horizontal direction. 

From a mathematical perspective, nutrient is repre¬ 
sented as a concentration field, the dynamics of which 
are governed by the the reaction-diffusion equation 


^^=V-(Ds(x)-V5(x))+rs(x), (1) 

where 5'(x) is the space (x) -dependent nutrient concen¬ 
tration, Dsi'x.) is the diffusion coefficient of the nutrient, 
and r 5 '(x) is the consumption rate of the nutrient by the 
bacteria. The rate of nutrient consumption, r 5 '(x), is 
related to the growth rate of the bacteria, dX/dt, via 


, , dS 1 dX 


( 2 ) 


where X(x) is the local biomass density, and Y^/s is a 
yield coefficient that describes the amount of nutrient 
required to produce one unit of biomass X. 

The growth rate of each cell is governed by the well- 
known Monod function 

^ 

ks ^ S ’ 

where jUrnax is the maximum specific growth rate of the 
bacteria, and ks is the concentration of nutrient, S, at 
which the growth rate is half maximal. The growth pa¬ 
rameters used in our simulations were taken from em¬ 
pirical and simulation studies on Pseudomonas aerugi¬ 
nosa^ assuming glucose to be the rate-limiting nutrient 
(see Table I). However our results are not sensitive to 
the detailed parameter choice, or to the choice of rate- 
limiting nutrient. Note that the growth rate parameters 
^x/Ni k'max^ and ks are the same for both the aggregate 
cells and the competitor cells (see Table I). 

From a practical point of view, in idynoMiCs the nu¬ 
trient concentration fields are assumed to be in pseudo 
steady-state with respect to biomass growth and there¬ 
fore the time dependence is removed from Equation 1 


0 = V.(D5(x).V5(x))-Fr5(x). (4) 

This equation is solved numerically for every global up¬ 
date of the bacterial population. In our simulations, we 
used a bulk nutrient concentration of 5.4 x 10“^ gL“^ 
(see Section S6) comparable with previous work [32-34]. 
Each of our simulations was run for a total time of 480 h. 
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Symbol 

Description 

Value 

Notes/ref 

[S] bulk 

Bulk concentration of limiting nutrient 

5.4 X gL“^ 

Within range of values from [32-34] 

Ux/AT 

Yield coefficient for Monod equation (Equation 3) 0.44 

Within range of values from [32, 35, 36] 

l~^max 

Maximum specific growth rate 

0.35 h-i 

Within range of values from [35, 37-39] 

Kn 

Half saturation concentration of nutrient 

3 X 

Within range of values from [20, 21, 33, 35-38, 40] 

7 

Density of biomass 

200 gL-i 

[20, 32, 35] 

Dg 

Diffusivity of glucose in water 

5.8“^ m^day“^ 

[41] 

Lx 

Dimension of system in horizontal direction 

1032 jim 

Ensures aggregates do not interact periodically 

Ly 

Dimension of system in vertical direction 

1032 jim 

Corresponds to the horizontal length 

Ldhi 

Thickness of diffusion boundary layer 

80 /xm 

Within range of values from [20, 21, 42] 


TABLE L Input parameters for biofilm simulations. 


in order to explore the long-term growth dynamics (see 
results and Section S7). 

Our complete parameter set is listed in Table L Using 
the parameters in Table I, our simulations produce spa¬ 
tially structured biofilms 200-300 /im in height after a 
simulation time of 480 h (see Figure 5, and Section S6). 


RESULTS 

Initial Aggregate Shape Determines Growth 
Dynamics 

To assess the growth dynamics of pre-formed aggre¬ 
gates in a biofilm, we tracked the number of progeny cells, 
A, produced by aggregates of different shape (Figure 
1(b)). We investigated three different aggregate shapes, 
characterised by the angle 0 (See Methods). The three 
angles, 0 = 5°, 90°, 180°, describe pre-formed aggregates 
that are initially arranged on the surface in either a 
spread, intermediately spread, or rounded manner. To 
investigate the effect of the surrounding unaggregated 
competitor cells (red cells in Figure 1(a)), we varied the 
density, p, of these cells between two extremes regimes of 
competition reported in this study: p = 0 cell pm“^, no 
competition; and p = 0.5 cell pm“^, high competition. 

Not surprisingly, aggregates grow better in the absence 
of surrounding cells on the surface regardless of their 
initial shape, i.e., for p = 0 cell pm“^. In this “non¬ 
competitive regime”, the initially spread aggregate pro¬ 
duces more progeny than the more rounded aggregate. 
This is evident in Figure 2, which shows representative 
initial aggregate configurations and the structures of the 
biofilms which they form after 480 hours. 

Competition from unaggregated competitor cells on 
the surfaces leads to more complex behaviour. Figure 
1(b) shows that, in the presence of strong competition, 
the spread aggregate produces more progeny over short 
times than the rounded aggregate. However, over longer 
times, the size of the population arising from the rounded 
aggregate is larger. 


For an aggregate in the presence of competing non¬ 
aggregating cells. Figure 1(b) points to two strategies 
for maximising progeny. For long-lived biofilms, progeny 
can be maximised by the aggregate adopting a rounded 
configuration, whereas if the biofilm is short-lived then 
it may instead be optimal for the aggregate to spread on 
the surface. 


Initial Aggregate Shape Affects Long-Term Biofilm 
Structure 

In our simulations the initial shape of the aggregate 
influences the long-term structure of the biofilm. Fig¬ 
ure 2 shows typical biofilm structures formed after 480 
h of growth, starting from aggregates that were initially 
spread on the surface (6> = 5°, left), rounded {0 = 180°, 
right) or partially spread {0 = 90°, centre), in the absence 
of competition from surrounding cells. It is clear that the 
spread aggregate covers much more of the surface during 
growth than its more rounded counterparts. 

In the presence of competition (Figure 3), we observe 
a marked difference in the structure of the biofilms that 
originate from spread aggregates (left panels) and from 
rounded aggregates (right panels). For the spread aggre¬ 
gate (a and c), the green section of biofilm that originates 
from the aggregate is structurally indistinguishable to 
that of the surrounding red biofilm that originated from 
the competing, unaggregated cells. In contrast, for the 
rounded aggregate (b and d), cells originating from the 
aggregate form a distinct “clump”, which is taller than 
the surrounding biofilm. When the density of competing 
(red) cells is high (Figure 3(d)), there is a cell-free gap 
around the growing clump that appears to be a result of 
nutrient depletion. 

It is clear from the biofilm structures shown in Figures 
2 and 3 that, even at very long times, the spatial structure 
of a biofilm can be affected by the initial spatial configu¬ 
ration of its founder cells (see also Section S6). While it 
might seem remarkable that apparently small changes in 
initial configuration can have dramatic effects on biofilm 
structure even after many cell generations, this effect is 
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FIG. 2. Initial aggregate arrangement affects biofilm morphology. Simulation snapshots of three bacterial aggregates initially 
arranged on the surface and the biofilms they form after 480 h: (a) Spread, 0 h. A zoomed in image is also shown to make the 
shape of the aggregate easier to resolve; (b) Semi-spread, 0 h; (c) Rounded, 0 h; (d) Spread, 480 h; (e) Semi-spread, 480 h; (f) 
Rounded, 480 h. 
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FIG. 3. Aggregate shape and neighbouring strain density 
affect biofilm morphology. Simulation snapshots of biofilms 
seeded from spread and rounded aggregates after 480 h growth 
in the presence of a low and high density inoculum of the 
competing strain: (a) ^ = 5°,p = 0.01 cell /xm”^; (b) 0 = 
180°,p = 0.01 cell pm“^; (c) ^ = 5°,p = 0.5 cell pm“^; (d) 
0 = 180°,p = 0.5 cell pm“^. 


in fact well-known in a different context. For initially 
flat biofilms, Dockery and Klapper showed theoretically 
that small inhomogeneities in initial configuration may be 
magnified into large “fingers” over the course of biofilm 
development [30]. This phenomenon is shown as a finger¬ 
ing instability and arises from the fact that an emerging 
protrusion (or finger) is elevated above, and thus depletes 
nutrients from the surrounding biofilm. This leads to 
positive feedback, in which the enhanced growth of the 
cells at the top of the instability is to the detriment of 
the surrounding cells below [43]. 

While Dockery and Klapper assumed that struc¬ 
tural inhomogeneities would arise spontaneously during 
biofilm growth, in our simulations such inhomogeneities 


are effectively created by the presence of the initial ag¬ 
gregates. The introduction of the rounded aggregate 
amongst the lawn of unaggregated cells on the surface 
at high competition leads to an instability in the biofilm 
structure that propagates as the biofilm develops. 


Nutrient gradients are important determinants of 
aggregate fate 



FIG. 4. Gells on the outside of the aggregates grow faster 
because they have greater access to nutrients, (a) Gell growth 
rate (p) distribution of the biofilm formed from the semi¬ 
spread aggregate in the absence of competition after 4h. (b) 
Gorresponding nutrient concentration field, [S]. 
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FIG. 5. Gradients in individual cell growth rates emerge in our simulated biofilms during growth. Gell growth rate distributions 
for the spread and rounded aggregates at after 480 h of growth: (a) 0 — ^ p — 0.0 cell (b) 0 — ^ p — 0.01 cell /xm”^; 

(c) ^ = 5°,p = 0.5 cell (d) 0 = 180°,/? = 0.0 cell (e) 0 = 180°,p = 0.01 cell (f) 0 = 180°,p = 0.5 cell 

pm“^. These distribution correspond to the conhgurations in Figures 2 and 3. Note that the gradient in cell growth rate is so 
large that a log scale is used for visualisation purposes. The green dashed lines represents an approximate boundary between 
the aggregate cells and the surrounding competing strain. 


It is well known that growth rate heterogeneities, re¬ 
sulting from nutrient concentration gradients, emerge 
during biofilm growth [44]. With this in mind, we tracked 
the growth rates of individual cells as a function of their 
position within the growing biofilm. Even in the very 
early stages of biofilm growth, we see heterogeneity in 
growth rates which emerge from (and influence) spatial 
gradients in nutrient concentration. Figure 4 illustrates 
this for a semi-spread aggregate {0 = 90^), after 4 h 
of growth, in the absence of competition. As expected, 
the cell growth rate is highly heterogeneous across the 
biofilm. Figure 4(a), with cells on the outside growing 
faster than those on the inside because they have better 
access to nutrients (Figure 4(b)). 

The growth rate heterogeneities shown in Figure 4(a) 
are amplified in the later stages of biofilm growth. Figure 
5 shows the spatial distribution of cell growth rates for 
biofilms arising from spread and rounded aggregates after 
480 h, in the presence and absence of competitor cells, 
for the same simulations as shown in Figures 2 and 3. 

In all cases we observe, as in previous work [45], two 
distinct regions of growth activity within the develop¬ 
ing biofilms: an outer layer of metabolically active cells 
and an interior region of inactive cells. These distinct 
regions arise because consumption by cells in the outer 
layer deprives cells in the inner layer of nutrients [45]. 
We also observe a large gradient in individual cell growth 
rate within the growing layer itself (note the logarithmic 
scale in Figure 5). The dynamics of the metabolically 
active layer determine the overall growth behaviour and 


structure of the biofilm. In Section S9 we show that the 
active layer of the rounded aggregate, unlike the spread 
aggregate, continues to expand in the presence of compe¬ 
tition; explaining why its total population becomes larger 
than that of the spread at longer times in Figure 1(b), 
and why it’s structures tend to fan outwards at the top 
(Figure 5(f)). 

Although it is well documented that nutrient gradi¬ 
ents arising during biofilm growth play an essential role 
in biofilm formation [8, 23, 29], so far few studies have 
investigated the effect of pre-formed bacterial clumps in 
this process; in particular how the initial arrangement of 
cells within a clump affects biofilm structure and develop¬ 
ment. Figures 4 and 5 show that the initial arrangement 
of cells on the surface can determine the shape and struc¬ 
ture of a growing biofilm because small initial differences 
in nutrient gradients become amplified as the biofilm de¬ 
velops. 


Competition for Nutrient Favours Rounded 
Aggregates 

Next we investigated how the fate of an aggregate, as 
measured by the average number of progeny of one of 
its cells, varies with aggregate shape. To this end, we 
computed the number of progeny cells, N, arising from 
the aggregate after a period of biofilm growth, relative 
to the initial number of cells in the aggregate, A'o, for a 
range of aggregate shapes, determined by 6>, at varying 
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FIG. 6. Average number of progeny, N/Nq, of aggregates 
defined by their surface-aggregate angle 0, the functional from 
of which changes with increasing density of competitor cells: 
(a) p = 0 pm celP^; (b) aggregate-medium interface length, 
s, as a function of 0; (c) p — 0.145 pm celP^; p = 0.5 pm 
cell“^. Vertical bars in represent the standard deviation from 
20 data points. 
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FIG. 7. Relative fitness as measured by N/Nq of rounded 
aggregates increases with competition. Rounded aggregates 
become favourable relative to spread aggregates with increas¬ 
ing density of competitor cells. 


levels of competition. For this analysis, we carried out 
long simulations (480 h of biofilm development), so that 
the ratio N/Nq reflects the long-time fate of the progeny 
of cells within the aggregate (see Section S7). 

Figure 6(a) shows N/Nq plotted as a function of the 
aggregate-surface angle, in the absence of competition 
from surrounding unaggregated cells. It is clear that the 
spread aggregate produces more progeny on average than 
the rounded one. 

In the previous section we saw that cells on the out¬ 
side of the aggregate have more access to nutrients even 
in the very early stages of biofilm growth. Thus, we might 
hypothesise that the growth advantage of the spread ag¬ 
gregate, in the absence of competition, is related to its 
larger surface area. Indeed, Figure 6(b) shows that N/Nq 
correlates closely with the interfacial area (or arc length) 
of the initial aggregate, s{0). We therefore conclude that 
the spread aggregate produces more progeny than the 
rounded one because the former has a greater surface 
area with the surrounding medium, providing greater ex¬ 
posure to nutrient in the initial stages of growth. The 
difference in initial structure between the spread and 
rounded aggregates therefore translates into significant 
differences in cell fate, even after many generations of 
biofilm growth. 

As the density of competition of unaggregated cells on 
the surface increases, however, a very different scenario 
emerges. Figures 6(c) and (d) show that cells in the 
rounded aggregate (large 0) produce more progeny, on 
average, than those in the spread aggregate. This is more 
evident in Figure 7, which shows N/Nq for the spread 
and rounded aggregates as a function of the density of 
competitor cells. This effect is also observed at higher 
nutrient concentration (see Section S8). 


Why does competition from surrounding unaggregated 
cells favour rounded over spread aggregates? Close in¬ 
spection of Figure 6 shows that, while the number of 
progeny produced by the spread aggregate decreases with 
increasing competition from surrounding cells (panels 
(c), and (d)), the number of progeny produced by the 
rounded aggregate remains rather constant. 

This finding can be understood by investigating how 
the fate of an individual cell within an aggregate depends 
on its initial spatial location. To this end, we tracked the 
number of progeny of each individual founder bacterium, 
as a function of its initial position within an aggregate. 
This constitutes a local, spatially-resolved version of the 
“fitness measure” N/Nq. Averaging our results over 20 
repeated simulations allowed us to generate a map show¬ 
ing the average number of progeny produced by indi¬ 
vidual cells within an aggregate, for initially spread and 
rounded aggregates. Figure 8. 

Figure 8 shows that the initial position within an ag¬ 
gregate indeed has a strong effect on cell fate. In the 
absence of competition from surrounding unaggregated 
cells, the most successful cells in the spread aggregate 
are those at the horizontal extreme edges; in the interior 
region of the aggregate, cell fate is more uniform (Figure 
8(a)). It seems likely that in this case, cells at the hori¬ 
zontal edges have an advantage because their progeny can 
expand in the horizontal direction, whereas the progeny 
of cells in the interior of the aggregate must compete 
with their neighbours within the aggregate for nutrients 
and space. The proliferation of the cells at the edges of 
the aggregate drives the lateral expansion of the growing 
biofilm which we observe in Figures 2(a) and (d). In con¬ 
trast, for the rounded aggregate (Figure 8(b)), cell fate 
is overall more heterogeneous within the aggregate, with 
the most successful cells being located around the Out- 


Spread 
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side surface of the aggregate. For the rounded aggregate, 
it appears that height is a relevant factor as well as the 
proximity to the aggregate surface. 

Figure 8(c) shows that, in the presence of competition, 
cells at the horizontal edges of the spread aggregate actu¬ 
ally do less well than those in the centre. The decreased 
fitness of these cells explains the inhibited lateral expan¬ 
sion observed in Figures 3(a) and (c). For the rounded 
aggregate, the most successful cells in the absence of com¬ 
petition are those at the top of the aggregate. In the 
presence of competition (Figure 8(d)), these cells, which 
are now highly localised at the top, are elevated above 
the level of the competitor cells and therefore are little 
affected by the increased competition for nutrients. The 
“fitness” cost associated with its smaller surface inter¬ 
face is compensated in the presence of competition by its 
height, since its top cells remain unchallenged by com¬ 
petitors with respect to nutrient access. 


DISCUSSION 

Given the tendency of many bacteria to aggregate, and 
the frequent observation of aggregates in diverse environ¬ 
mental situations [5, 46, 47], it seems likely that natural 
biofilms are often initiated from pre-formed aggregates. 
Despite this, the role of pre-formed aggregates in biofilm 
development has, to our knowledge, not yet been ad¬ 
dressed. In this paper, we have investigated the fate 
of pre-aggregated cells during biofilm formation, using 
individual-based simulations. Our study shows that an 
initial aggregate can have a significant and long-lasting 
effect on biofilm spatial structure, even after many gen¬ 
erations of cell growth. Focusing on the role of aggre¬ 
gate shape, we find that, in the absence of competition 
for nutrients from surrounding cells, an aggregate that is 
initially spread on the surface is favoured over one that is 
initially rounded even over long periods of biofilm devel¬ 
opment. This is likely to be because the spread aggregate 
has initially a larger surface over which it can absorb nu¬ 
trients, giving it an initial growth advantage that is then 
maintained as the biofilm grows. 

Strikingly though, our results change qualitatively in 
the presence of competition from surrounding, unaggre¬ 
gated cells. When this competition is strong, although 
the spread aggregates still grow faster in the early stages 
of biofilm development, rounded aggregates become more 
successful (produce more progeny) as the biofilm devel¬ 
ops over longer times. This effect appears to arise from a 
trade-off between height (as nutrients diffuse from above) 
and exposed surface area. In the absence of competi¬ 
tion, surface area is more important than height, and the 
spread aggregate is favoured. However, in the presence 
of competition, height becomes more important, since 
cells at the top of the aggregate can avoid competing for 
nutrients with the surrounding competitors. Since the 
rounded aggregate is taller than the spread aggregate, 
it gains a “fitness” advantage under conditions of strong 


competition that is only realised after long times. 

Bacterial biofilm formation is a complex phenomenon 
which involves a plethora of biological mechanisms in¬ 
cluding cell motility [48], EPS production [49, 50], 
metabolic and other phenotypic differentiation [9, 47, 51], 
and cell-cell interactions such as quorum sensing [52-54]. 
In our simulations, almost all of this biological complex¬ 
ity has been neglected; our model takes account only 
of nutrient gradients established by cell consumption, 
nutrient-dependent growth, and competition among cells 
for space. Nevertheless this simplistic approach produces 
biologically interesting, and potentially testable, predic¬ 
tions. In particular our simulations predict that being 
initially spread on a surface is a better strategy for a 
bacterial aggregate in the absence, but not in the pres¬ 
ence, of competition. Understanding how further bio¬ 
logical complexity might affect this picture would be a 
very interesting topic for further work. Another avenue 
worth investigating would be the effects of biofilm erosion 
and the subsequent detachment of cells. Here, our simu¬ 
lations have not included the effects of fluid flow, which 
among other effects, may flatten the biofilm by detaching 
protruding cells. 

How might aggregates of different shape arise in na¬ 
ture? It is well known that bacterial interactions with 
surfaces can vary greatly depending both on the physi¬ 
cal and chemical properties of the surface [55-57], and 
on bacterial phenotypes such as EPS production and the 
presence of surface appendages. It is therefore very likely 
that aggregates formed from bacteria of different taxa or 
strains, landing on different surfaces, might adopt dif¬ 
ferent configurations. Eor example, certain bacteria pro¬ 
duce surfactant which can alter the morphology of a de¬ 
veloping biofilm and allow them to expand over surfaces 
more efficiently [58, 59]. 

Our work has been inspired by the observation that 
bacterial aggregates often form in the planktonic phase 
[11-13]. Aggregates are known also to form via the 
detachment of bacterial clumps from a mother biofilm 
[15, 60, 61]; should such aggregates land on a pristine 
surface, similar phenomena to those discussed here would 
be expected to arise. Moreover, our results could also be 
relevant to aggregates that form on the surface itself. 
In the classical picture of P. aeruginosa biofilm devel¬ 
opment, individual cells land on a surface, upon which 
they migrate and proliferate to form small aggregates 
(i.e. microcolonies). Surface-induced motility mecha¬ 
nisms [62, 63] such as twitching [64], crawling [65] and 
walking [66] have been implicated in this process. Once 
formed, such small aggregates would compete with sur¬ 
rounding cells for nutrients in much the same way as the 
pre-formed aggregates that we investigate in this paper. 

This study has focused on a single pre-formed aggre¬ 
gate seeding the surface and competing with initially un¬ 
aggregated cells during biofilm formation. To further un¬ 
derstand biofilms in nature, this work should be extended 
to investigate competition between multiple aggregates 
arranged on the surface, and competition between ag- 
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FIG. 8. Distribution of fittest cells varies with aggregate shape. 2D histograms representing the number of progeny, N^ produced 
(480 h) by individual bacteria as a function of their initial location in the spread and rounded aggregates in the absence and 
presence of competition: (a) ^ = 5°,p = 0.0 cell /xm”^; (b) 6 — 180°, p = 0.0 cell pm“^; (c) ^ = 5°,p = 0.5 cell pm”^; (d) 
0 — 180°,p = 0.5 cell pm“^. Note that these distributions were averaged over 20 trajectories for each aggregate. Note that the 
gradient in the number of progeny is so large that a log scale is used for visualisation purposes. 


gregates and mixed strains of bacteria, i.e., strains with 
different growth rates. 

Recently cooperation within clumps of aggregated cells 
has been suggested to be a major stepping stone in the 
evolution of multicellularity [26, 27]. Our study thus also 
hints that interesting social interactions might arise be¬ 
tween cells within an aggregate. For all aggregate shapes, 
we observe heterogeneity in fitness among cells within 
the aggregate. This is particularly pronounced for the 
rounded aggregate, where cells at the top are strongly 
favoured while those in the centre of the aggregate hardly 
proliferate. Based on arguments recently put forward 
by West and Biernaskie [26, 27], one might predict that 
rounded aggregates would be favourable under condi¬ 
tions where cells within the aggregate are closely related, 
whereas spread aggregates, in which fitness differences 
between cells are less pronounced, might form where cells 
are less closely related. This leads to interesting fur¬ 
ther questions, e.g., when a rounded aggregate initiates 
biofilm growth, does the majority of cells in the aggre¬ 
gate “sacrifice” their future progeny in favour of their kin 
at the top? This idea supports previous suggestions that 
height plays a crucial role in competition within biofilms 
[32]. While previous work pointed to EPS production 
as a means to push progeny cells above the surrounding 
competitors [32, 67], our work shows that aggregate for¬ 
mation also provides a means to this end. Such a picture 
raises new questions about the evolutionary implications 


of bacterial aggregation. 
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Supplementary Material 


Shaping the Growth Behaviour of Biohlms 
Initiated from Bacterial Aggregates 


SI Generating 2D bacterial aggregates 

To generate bacterial aggregates of the appropriate shape, circular segments 
were extracted from simulation conhgurations of pre-grown biohlms as illus¬ 
trated in Figures SI and S2. Figure SI (a) shows a circle (red) with radius R 
and segmental region dehned by the white area between the black horizontal 
line and the perimeter of the circle. The “surface-aggregate angle”, 0, dehnes 
the angle between the black horizontal line (surface) and the tangent line of 
the circle at the point of contact with the surface. 

The area, A, of a circular segment of radius R and with surface-aggregate 
angle of 0, is given by [1] 


^ = l(26»-sin2^)i?2. (SI) 

The corresponding arc length, 5, of an aggregate dehnes the initial interface 
between the bacterial aggregate and the surrounding nutrient medium, and 
is computed according to 

s = 2Re. (S2) 

To ensure our aggregates contained the same number of cells (~100), the 
area in Equation SI was hxed at 1257 //m^. 

The procedure for generating the starting aggregate conhgurations was 
performed as follows: 


1 



1. An in silico biofilm was grown for a simulation time of 7 days. This 
simulation was initiated with 200 cells randomly distributed on the 
surface and a nutrient concentration of 5.4x10“^ g was used. A 
typical biohlm resulting from such a simulation was approximately 800 
//m in height with some finger-like projections such as those in Figure 
S2. 

2. A point, P, was defined in the middle of the biofilm (Figure SI and 
S2). The point P is the reference point from which the origin, O, of 
a circle is computed. The point P was selected such that the density 
of cells within any segmental area was relatively uniform. Bacterial 

aggregates of shapes 0 = 5,10,15, 20, 25,., 180^ with fixed area of 

1257 //m^ were then generated as in Figure S2. 

3. For each 0, the radius of the corresponding circle was computed by 
solving Equation SI for R. Using the radius, P, of this new circle, we 
computed the distance, X, to move the origin of this new circle above 
or below the point P according to 

X = R sina, (S3) 

where 

a = 90-9. (S4) 

4. With the new origin and radius defined, we then extracted all bacteria 
within the radial distance from O. Note those bacteria located below 
the horizontal line (blue region) were excluded. The corresponding 
segment gave rise to a bacterial aggregate with a shape determined by 
the value of 0. 

5. The vertical component of the point P was then subtracted from all 
bacterial coordinates so that the base of the aggregate was located on 
the surface as shown in Figure 1(a) of the main manuscript. 

S2 Generating configurations of competing “sin¬ 
gle cells” in 2D 

The width of the surface available to the surrounding competitor cells de¬ 
pends on the surface coverage of the aggregate. As surface coverage of the 
aggregate increases, 0^0, the amount of available surface either side of the 
aggregate decreases. Surrounding cells were inserted at random positions 
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Supplementary Figure SI: Schematic representation of the geometry of ag¬ 
gregates in (a) 2D and (b) 3D. The light blue region represents the area not 
included when computing cells within a distance R from the origin, O. The 
angle 6 is the contact angle of the aggregate with the surface. 



Supplementary Figure S2: Generating aggregates with different geometries 
but containing the same number of cells. Using pre-grown 2D biofilm con¬ 
figurations, green, aggregates of the desired shape can be generated using 
the geometric relationships from Figure Sl(a). (a) Schematic showing how a 
spread aggregate is generated, (b) Schematic showing how a more rounded 
aggregate is generated. 
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on the surface excluding the region occupied by the aggregate. To ensure a 
constant density of cells surrounding the aggregate, the number of cells in 
these regions were varied according to the width of surface available. 

53 Generating 3D aggregates 

To generate aggregates in our 3D simulations, we performed an analogous 
process to that used for the 2D case. Conhgurations of spherical caps were 
transplanted from pre-grown biohlms using the geometrical relationships il¬ 
lustrated in Figure Sl(b). The volume of a spherical cap is [2] 

Vcap = -7ri?^(2 — Ssina + sin^a), (S5) 

o 

and the corresponding surface area of the cap is given by 

S = 27cRh, (S6) 

where the height, /i, is given by 

h = R-X (S7) 

A spherical cap volume, Vcap^ of 5575 //m^ was used for all values of 0, giving 

approximately 118 bacteria in each aggregate. 

54 Generating configurations of competing “sin¬ 
gle cells” in 3D 

To surround the 3D aggregates with neighbours, we seeded the surface sur¬ 
rounding the aggregate with single cells. The number of these competing 
cells per square fim gave the density of cells on the surface. Cells that lay in 
the region dehned by the radius of the spherical cap, a, were removed from 
this conhguration. The aggregate conhguration and the surrounding cell 
conhguration were then combined to produce the initial state of the system. 

55 Multiple simulation runs in 2D and 3D 

For our 2D simulations, four different aggregate configurations were gener¬ 
ated for each value of 0, and for each of these configurations, 5 simulations of 
480 hours were performed using a different initial distribution of surrounding 
cells. This gave rise to a total of 20 simulations for each value of 6. 
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In 3D, only the spread aggregate, 0 = 5^, and the rounded aggregate, 6 = 
180^, were simulated. Three different initial configurations were generated 
for each value of 0, and for each of these configurations 3 simulations of 72 
hours were performed, each using different initial distributions of surrounding 
single cells. This gave rise to a total of 9 simulations for each value of 6. 

Note that, in both the 2D and 3D simulations, the competitor cells on 
the surface and the aggregated cells are identical in their growth parameters 
(Table 1, main manuscript). 


S6 Exploration of Biofilm Structure at Low 
and High Nutrient Concentration With¬ 
out Aggregates 

To investigate the effect of nutrient concentration in general in our simu¬ 
lations, we simulated 2D biofilms in the absence of initial aggregates, at 
several different nutrient concentrations and various densities of cells, which 
were initially located at random across the surface. Figures S3 (a) and (b) 
show the resulting biofilms formed after 120 h at a nutrient concentration of 
5.4 X 10“^g for seeder cell densities of 0.01 cell jivcT^ and 0.5 cell jivcT^. 
Clearly, the biofilms formed at this nutrient concentration are significantly 
thicker than those formed at lower nutrient concentration (Figures S3(c) to 
(f)). Comparing (a) and (b), we see that a low density of seeder cells leads 
to a less uniform biofilm morphology with finger-like projections. 

After 120 h, at low nutrient concentrations (5.4 x 10“^g L“^), the biofilms 
remain very thin (Figures S3(c) and (d)). At this low concentration, the 
initial number of seeder cells has a marked effect on the final structure of 
the biofilm. A high density of seeder cells yields a biofilm with a uniform 
morphology whereas a low density of seeder cells yields a biofilm with a 
spatial structure that is clearly dependent on the initial configuration of the 
seeder cells. Figures S3(e) and (f) show that this effect persists at long times. 

From Figures S3(a) and (e), it is clear that both longer times and/or 
higher nutrient concentration lead to larger biofilms. Selecting a bulk nutrient 
concentration and timescale to investigate the shape dependent fitness of 
aggregated clusters of bacteria is therefore non-trivial. In our simulations, 
we used the lower nutrient concentration of 5.4 x 10“^g (see Table 1 in 
the main manuscript), comparable with previous biofilm simulations [3, 4, 5]. 
We also chose to run the simulations for 480 h in order to explore the long 
term growth dynamics, and to generate populations that were large enough 
to obtain good statistics. Timescales of 480 h have also been used in previous 


5 




Supplementary Figure S3: Nutrient concentration and initial density of 
seeder cells affect biofilm morphology in the absence of initial aggregates: 
(a) Nutrient concentration = 5.4 x 10“^g Density of seeder cells = 0.01 
cell . (b) Nutrient concentration = 5.4 x 10“^g Density of seeder 
cells = 0.5 cell (c) Nutrient concentration = 5.4 x 10“^g Density 

of seeder cells = 0.01 cell (d) Nutrient concentration = 5.4 x 10“^g 

Density of seeder cells = 0.5 cell (e) Nutrient concentration 

= 5.4 X 10“^g Density of seeder cells = 0.01 cell //m“\ 480 h. (f) Nutri¬ 
ent concentration = 5.4 x 10“^g Density of seeder cells = 0.5 cell 
480h. 


simulation studies [3, 4]. After 480 h growth using a bulk concentration of 
5.4 X 10“^g our simulations produce biohlms that are approximately 
200-300 jim. in height (Figures S3(e) and (f)). 
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Supplementary Figure S4: Without competing cells on the surface, the func¬ 
tional form of the N/Nq curves remains relatively constant over time. Note, 
however that the scale does change because more cells are produced as time 
increases: (a) 56 h. (b) 120 h. (c) 240 h. (d) 480 h. Error bars represent the 
standard deviation from 20 simulations. 


S7 Exploration of Aggregate Fate at Differ¬ 
ent Times 

In the main manuscript, we computed N/No^ as a function of the aggregate 
surface angle, 0, after 480 h. Here, we show that qualitatively similar results 
are obtained regardless of the time at which htness is measured. Figures S4 
shows N/Nq as function of 6 in the absence of competing cells (no red cells 
in Figure 1(a) of the main manuscript) at 4 different times. Evidently, the 
functional form of the curves does not change significantly with time. 

In the highly competitive regime where we include competitor cells (red 
cells in Figure 1 of main manuscript) at a density of 0.5 cells the 

situation is more complex. Here, the functional form of the curves changes 
during the course of the simulation (Figure S5). For early times, the curves 
are similar to the case without competition, suggesting that shape-dependent 
nutrient access is the still dominant factor in the early stages. As time 
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Supplementary Figure S5: With competing cells on the surface at a density 
of 0.5 cell the functional form of the N/Nq curve changes over time: 

(a) 56 h. (b) 120 h. (c) 240 h. (d) 480 h. Error bars represent the standard 
deviation from 20 simulations. 


increases, the rounded aggregates become increasingly “htter” relative to 
the spread aggregates, until at 480 h, all aggregates with a surface aggregate 
angle greater than 75^ are htter than the spread aggregate. From this we 
can conclude that the initial height advantage of the rounded aggregate only 
becomes important at longer times. 


S8 Difference in Fate Between Spread and 
Rounded Aggregates at Higher Nutrient 
Concentration 

Here, we investigate the fate of the spread and round aggregate at a higher 
nutrient concentration, 1 x 10“^ gL“^. Figure S6 shows that our conclusions 
remain valid at this nutrient concentration: at low density of seeder cells there 
is no difference in fate between the spread and rounded aggregates, whereas 
at higher density of initial neighbours, the rounded aggregate is favoured 
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Supplementary Figure S6: At higher nutrient concentration, 1 x 10-2gL-i, 
rounded aggregates are more favourable for growth. At low density of seeder 
cells (0.01 cell there is no signihcant difference in N/Nq between the 

spread and rounded aggregates. At high density (0.5 cell the rounded 

aggregate is again produces more progeny per initial cells relative to that of 
the spread. Note that these data were generated from simulations represent¬ 
ing 480 h of biofilm growth. 


over its initially more spread counterpart. Note that these simulations were 
again run for 480 h. 


S9 The Effect of Aggregate Shape on the Growth 
of the Active Layer 

In the main manuscript we discussed cell growth rate heterogeneity within 
the biofilm and the presence of a well-defined active layer in the growing 
aggregate. The depth of the active later has been shown to be very important 
in determining biofilm structure and the degree of segregation amongst cell 
groups within the community [6]. It has also been discussed in extensive 
mathematical detail by Dockery and Klapper [7]. Here we look at the growth 
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Supplementary Figure S7: Dynamics of the growing layer determine aggre¬ 
gate growth behaviour. The number of growing cells in the spread and 
rounded aggregates at various competitive regimes. For clarity the error 
bars, representing the standard deviations, are only shown for the hnal data 
points. The standard deviations at these points are maximal. 


of the active layer. 

In the analysis that follows, we simply dehne the active layer as the 
region of the growing aggregate in which all cells grow with a growth rate 
greater than 1 fg h“^ (see Figure 5 of main manuscript), and we assess the 
number of cells in this layer. This choice is arbitrary; what is important is 
the qualitative difference in behaviour between the aggregate shapes. 

Figure S7(a) shows the number of active cells within the spread {6 = 
5^) and rounded {6 = 180^) aggregates at varying degrees of competition; 
dictated by the density of surrounding unaggregated cells on the surface. In 
the absence of competition, the number of active cells within the spread and 
rounded aggregates increases with time. Initially the number of active cells 
in the spread aggregate grows faster than that for its rounder counterpart, 
however for t > 100 h, the number of active cells in both aggregates increases 
with similar rates (slope of red curves). This linear increase determines the 
super-linear increase in the corresponding aggregate populations in Figure 
1(b) of the main manuscript: a linear increase in the number of active cells 
translates to a super-linear increase in the total aggregate population. 

In all competitive regimes, the number of cells in the growing layer of the 
spread aggregates becomes constant after ^100 h. A constant population of 
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active cells in the growing layer translates to linear growth behaviour of the 
total aggregate population (compare solid blue curve with the dashed red 
curve in Figure 1(b) of the main manuscript). For the rounded aggregates, 
this transition to the constant regime does not occur within the timescales 
of our simulations, and the layer continues to grow monotonically. 

Comparing the spread and rounded aggregates (solid curves with dashed 
curves), with competition, we see that growth of the active layer always in¬ 
creases faster for the spread aggregate before levelling off after 100 h. How¬ 
ever, later in the simulation, the number of active cells in the rounded ag¬ 
gregates becomes greater. The time at which this crossover occurs increases 
with increasing levels of competition. The increasing growth of the active 
layer in the rounded aggregate at high competition, p = 0.5 cell ex¬ 

plains why its number of progeny becomes larger than that of the spread 
aggregate in Figure 1(b), and why its structures tend to fan outwards at the 
top (Figure 4(d)). 


SIO Aggregate Simulations in 3D 

The results presented in the main manuscript are for 2D simulations, for 
reasons of computational feasibility. To validate our results further, we also 
ran some simulations in 3D. Figure S8 shows simulation snapshots after 72 
h for the spread and rounded aggregates in the presence of low and high 
densities of competitor cells. As in our 2D simulations (Figures 3(b) and 
(d) of the main manuscript), the rounded aggregate persists as a dominant 
structural feature at both low and high densities of competitor cells. The 
spread aggregate on the other hand becomes swamped by the surrounding 
competitors and after 72 h is structurally indistinguishable from the rest 
of the biohlm (as in our 2D simulations. Figures 3(b) and (d) of the main 
manuscript). 

N/No for the spread and rounded aggregates are plotted in Figure S9 and 
is analogous to the data shown in Figures 8 of the main manuscript for the 
2D simulations. As in the 2D case, the number of progeny per initial cell 
of the aggregates decreases with increasing density of surrounding cells. In 
the absence of competition, the spread aggregate again produces signihcantly 
more progeny than the rounded aggregate, however as the competition on the 
surface is increased, the number of progeny produced by cells in the rounded 
aggregate increases relative to that of the spread. 


11 




Supplementary Figure S8: 3D simulation snapshots showing initial and hnal 
conhgurations of the spread and rounded aggregates at low and high density 
of surrounding cells, (a) Spread aggregate, Density of seeder cells = 0.002 
cell (b) Spread aggregate, Density of seeder cells = 0.03 cell yum“^. (c) 

Rounded aggregate. Density of seeder cells = 0.002 cell (d) Rounded 

aggregate. Density of seeder cells = 0.03 cell 
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Supplementary Figure S9: N/Nq of the rounded aggregate with respect to 
the spread aggregate increases with the density of neighbouring cells on the 
surface. In 3D, the trend is similar to that observed in 2D. The spread 
aggregate is signihcantly favoured when there is no competition on the surface 
however with increasing density the rounded aggregate is more favourable for 
growth. 
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